Extended Fractons and Localized Phonons on Percolation Clusters 
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We investigate the localization behavior of vibrational modes of infinite percolation clusters above 
the critical concentration in two and three dimensions using the method of level statistics. We find 
that all eigenstates are localized in d — 2, including the low frequency phonon states. In d = 3 we 
obtain evidence for a localization-delocalization transition. But contrary to the common view this 
transition occurs for frequencies above the phonon fracton crossover giving rise to a new regime of 
extended fracton states. 

PACS numbers: 63.50,+x, 64.60.Ak, 63.20.Pw 
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Percolation represents one of the standard models for 
disordered systems. Its applications range from trans- 
port in composites and in amorphous and porous media 
to the properties of branched polymers and gels, for re- 
cent reviews see e.g. One intriguing aspect is the 
vibrational dynamics of percolation networks. Since the 
pioneering discovery of localized lattice vibrations ( "frac- 
tons") by Alexander and Orbach much work has 
been devoted to the study of their propertie s [pH l4[ and 
the investigation of various applications |l5|-|l8|]; see Jl9| 
for a recent review. 

It is well accepted, that all vibrational excitations are 
localized at the percolation threshold, and the density of 
states (DOS) scales as p(w) ~ w ds_1 , where d s « 4/3 is 
the spectral dimension j^-||,|2(|. Above criticahty, there 
exists a crossover frequency uj^: for frequencies larger 
than oj£, g{u) ~ u ds ~ 1 , while below o^, the DOS scales 
in the normal fashion, g(uj) ~ u) 1 . The meaning of 
the "fracton-phonon crossover" at uj^ and in particu- 
lar its implication on the localization behavior has been 
discussed in a controversial way. Analogies with ran- 
dom walks seem to indicate that the crossover charac- 
terizes a transition between localized fractons and ex- 
tended phonons. This view has been mostly accepted 
for three-dimensional systems ( [|i]||||,||,|l3|,|lG 19 1, but see 
also [|l5|), but not yet been proven. Analogies with the 
Anderson model , in contrast, suggest the absence of a 
localization-delocalization transition for two-dimensional 
systems. But since these analogies are not exact, even the 
value of the critical dimension has not been established 
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It is the purpose of this Letter to clarify the localiza- 
tion properties of the vibrational excitations. We em- 
ploy the method of level statistics which has been used 
successfully to determine the localization-delocalization 
transition (LDT) in the Anderson model [§2|-|^|. We 
find that the localization properties of the vibrations are 
not related to the fracton-phonon crossover (FPC) at uj^. 
Our results suggest that (i) in d = 2 all states are lo- 
calized and (ii) the LDT in d — 3 occurs at a critical 



frequency u> c well above u)£ . Accordingly the vibrational 
modes in d — 3 are extended phonons below u>£ and lo- 
calized fractons above tu c , while the regime between u)£ 
and lo c is characterized by extended excitations with a 
fracton- type dispersion ("extended fractons"). 

We consider the incipient infinite site percolation clus- 
ter on a square lattice and a simple cubic lattice, both 
above the critical concentrations p c ~ 0.59275 and ~ 
0.3116, respectively. On length scales below the correla- 
tion length £ ~ \p — p c | _I/ , the cluster is self similar and 
characterized by the fractal dimension df. Above £, the 
cluster appears homogeneous. The correlation exponent 
v, the fractal dimension df, and the spectral dimension d s 
depend on the space dimension d: v = 4/3, df — 91/48 
and d s ~ 1.3173 in d = 2 pg|; v ~ 0.875, d f ~ 2.524, 
and d s ~ 1.328 in d = 3 

We assume that equal masses M are placed on each oc- 
cupied site and nearest neighbor cluster sites n and n + S 
are coupled by equal (scalar) force constants k„ tn+ s- In 
this case, different components of displacements decou- 
ple and we obtain the same equation of motion for lattice 
vibrations for all components u n (t), 



, r d 2 

M - 2 u n [t) 



} J k n ,n+s [u n+ s(t) - u n (t)] 



(1) 



where the sum runs over all nearest neighbor sites n + 8 
of site n. The standard ansatz u n (t) = u nMJ exp(— iuti) 
leads to the corresponding time independent vibration 
equation. Setting AI — k n:n+ s = 1 (which determines 
the unit of the frequency) we obtain 



■ oj u r , 



(2) 



which represents an eigenvalue equation that can only be 
solved numerically. 

For obtaining an analytical estimate of the way the 
characteristic length in the vibration problem depends 
on the frequency lo and on the correlation length £, one 
can employ the analogy between the corresponding 
diffusion equation and the vibration equation, Eq. (pi). 
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In the diffusion problem, the characteristic length is the 
root- mean-square displacement R(t) — {r 2 {t)) x / 2 of a 
diffusing particle (random walker). At criticality, R{t) 
scales as R(t) ~ t 1 /^ with the random walk dimension 
d w = 2df/d s . Above p c , the finite correlation length £ 
gives rise to a finite 'correlation' time t% ~ £ dw . Below t$, 
the random walker discovers the fractal part of the infi- 
nite cluster and R(t) ~ i 1 /^, while above i^, on length 
scales where the cluster is homogeneous, diffusion is nor- 
mal and R(t) ~ i 1 ^ 2 . Since the (discretized) diffusion 
equation reduces to the vibration equation, Eq. (|I|), when 
the first time derivative is substituted by a second time 
derivative, we can obtain a corresponding characteristic 
length A(w) for the vibration problem by substituting t~ l 
in R(t) by {l/tf = u? . This yields |,| 



A( w ) 



for uj ^> ■ - 
for a; -C Wf- 



-1/2 



t-dw/2 



(3) 



The expression for the crossover frequency u)£ that sepa- 
rates the fracton from the phonon regime can be written 
as 



uj^ =u (p- p c ) 



d v u/2 



(4) 



where the proportionality factor loq depends on the de- 
tails of the underlying lattice. 

The scaling for the vibrational DOS g(u>) can be de- 
rived in an analogous way. By construction, the crossover 
frequency separating the fracton regime with g(uj) ~ 
w d s -i f rom ^ n e phonon regime with g(u>) ~ Lo d ~ 1 is ex- 
pected to be identical with uj^ from Eq. (0). We have 
determined luq from the DOS for percolation clusters on 
the square lattice and the simple cubic lattice for several 
concentrations above p c . The crossover from the fracton 
to the phonon regime is smooth (see also flj|). We find 
u>q « 15 in d = 2 and loq ks 6 in d = 3 for the parameter 
in Eq. (^). Both numbers are in good agreement with 
earlier (less extensive) numerical calculations of the DOS 

HI- 

In order to see the significance of the characteristic 
length A described by Eq. (J3J) , we have numerically solved 
Eq. (0) using the Lanczos algorithm. Figure [l] shows the 
wavelength X(u>) that we calculated by Fourier transform 
of the eigenstates at several frequencies and concentra- 
tions in d — 2 and in d — 3. The data nicely confirm that 
the wavelength A is the characteristic length A discussed 
above. At the FPC, the wavelength crosses over from 
fracton behavior for high frequencies to phonon behavior 
for low frequencies according to Eq. (^). The values of 
the crossover frequency u>£ are well described by Eq. (^). 
In addition, the prefactor loq is exactly the same as for the 
crossover frequencies that were previously determined for 
the vibrational DOS g(u>), since the crossover occurs ex- 
actly at ui = uj^ in Fig. |l|. Thus, the characteristic length 
A can be identified with the wavelength of the modes. 
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FIG. 1.: Scaled wavelength \(uj)(p — p c ) u versus scaled fre- 
quency u/ui£ for site percolation clusters in d — 2 and d = 3 
with U£ from Eq. (^) and uq determined from the scaling 
behavior of the DOS. The upper curve (d = 3) has been 
shifted by a factor of three. The symbols correspond to the 
concentrations of occupied sites, for d = 2: p = 0.61 (□), 
0.63 (O), 0.66 (A), 0.70 (v), 0.75 (+), 0.82 (*), and 0.90 
(O); for d = 3: p = 0.32 (□), 0.35 (O), 0.37 (A), 0.40 (v), 
0.44 (O), and 0.50 (*). The straight lines are fits to the data 
with slope —1 in the phonon regime (uj < oj^) and — 2/d w in 
the fracton regime (uj > uj%) according to Eq. (|^). In d = 3, 
we used the effective value d w = 3.4 that governs diffusion 
on length scales of the order of the correlation length con- 
sidered here. The asymptotic value d w = 3.8 applies only to 
considerably larger length scales. The figure shows that uj£ 
is the FPC frequency for the wavelength. 

Next we address the localization behavior of the vi- 
brational excitations. We want to know if there is an 
intimate relationship between the FPC and the localiza- 
tion behavior, such that, for example, phonons are ex- 
tended and fractons are localized, as one might infer from 
the analogies with random walks discussed above. This 
analogy suggests a LDT both in d = 2 and d = 3. On 
the other hand there is no LDT in the somewhat related 
Anderson model in d = 2 and computer simulations for 
fractons and electrons (in tight-binding approximation) 
on percolation clusters at p c have shown that the local- 
ization behavior of fractons and electrons is very similar 
for corresponding frequencies and energies [0. While 
in the quantum mechanical description the localization 
of the electrons is caused by scattering and interference, 
the vibrational modes are waves already in the classical 
description and may become localized by scattering and 
interference like the electrons. 

In the following we address the questions: (i) Are the 
phonon states in d = 2 (for uj < uj^) extended or localized 
states? (m) Does the localization-delocalization transi- 
tion (LDT) for vibrations in d = 3 occur at the same 
frequency as uj^7 For answering these questions, we have 
applied the method of level statistics, which proved to be 
a powerful tool for the electronic case recently |22]-|24|] . 
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In conducting disordered electronic systems, the energy 
spacing distribution P(s) of consecutive eigenvalues (lev- 
els) Ei shows the universal random matrix theory re- 
sult [p5 26 1 , which is well approximated by the appropri- 



ate Wigner surmise, P(s) = (tt/2)s exp(— us 2 /A). Here, 
s = \Ei + i — Ei\/A where A is the mean level spacing in 
the energy interval considered. For localized states the 
uncorrelated eigenvalues are described by the Poisson dis- 
tribution, P(s) = cxp(— s). In contrast to P(s = 0) = 1 
for localized states, the probability that two eigenval- 
ues are close in energy decreases to zero for extended 
states, because of the level repulsion due to the overlap 
of the corresponding eigenstates. For finite systems, the 
shape of P(s) is in between the two limiting cases and ap- 
proaches one of them with increasing system size. P(s) is 
system size independent only directly at the transition. 
Hence, the knowledge of the nearest neighbor spacing 
distribution P(s) for several system sizes tells whether 
the eigenstates are localized or extended in the limit of 
infinite system size. 
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FIG. 2.: Second moments 7o(w) of the level spacing distri- 
bution P(s) for vibrational modes of 3d percolation clusters 
with site concentration p = 0.37. The symbols denote dif- 
ferent system sizes L (radius of the clusters in topological 
space), L = 50 (□), L — 62 (O), L — 75 (A), L — 90 
(V), L = 110 (O), and L = 150 (+), corresponding to cu- 
bic systems with linear sizes ranging from 30 to 107. The 
correlation length was determined to be £ = 3.3. The LDT 
occurs for u c = 0.24. 

Since it is tedious to compare the P(s) distributions 
for several system sizes and other model parameters, we 
calculate the size dependence of the quantity 1$ = (s 2 )/2, 
which is related to the second moment of the level spac- 
ing distribution (s 2 ) = J °° s 2 P(s) ds. If the eigenstates 
are localized, Iq increases monotonously with increasing 
system size approaching the Poisson limit I Q ° C = 1 for 
infinite system size. If the modes are extended, Iq de- 
creases monotonously until it reaches the random matrix 
theory limit 7g xt = 0.637. At the transition I™' is scale 
invariant, but its precise value depends on the bound- 



ary conditions. Other moments of the distribution give 
similar results. 

We apply the method of level statistics to vibrations 
of incipient infinite percolation clusters at several con- 
centrations ranging from p = 0.61 to p — 0.8 in d = 2, 
and from p = 0.32 to p = 0.5 in d — 3. Up to six 
different cluster sizes are compared for each concentra- 
tion. The clusters were generated by the Leath method 
and the eigenvalues were computed with a Lanczos al- 
gorithm. For the vibrational modes, the level spacing is 
s = — uj 2 \/A since the eigenvalue uj 2 in Eq. (||) cor- 
responds to the energy eigenvalue E for electrons. Rep- 
resentative results are shown in Fig. |^ and in the inset of 
Fig. | 

In Fig. U the LDT can clearly be observed for vibra- 
tions of clusters in d = 3 at p = 0.37. For low frequencies, 
u> < uj c , Iq(u!,L) decreases with increasing system size L, 
indicating extended modes, while it increases for lu > lo c , 
as expected for localized modes. Directly at the critical 
frequency uo c « 0.24, Iq(u,L) is independent of L. This 
value of the LDT frequency lo c is well above the corre- 
sponding FPC frequency w 0.053 that we obtained 
from Eq. @) with wo ~ 6. Consequently, we discover a 
novel intermediate regime of extended fracton modes for 

LU£ < LU < LU C . 

Using plots similar to Fig. || we determined the LDT 
frequency co c (p) for nine concentrations p. All values of 
uj c (p) are well above the respective FPC at oj^(p). Our 
results for d — 3 are summarized in Fig. ^ which shows 
the phase diagram. The critical frequencies uj c (p) can be 
fitted by w c (p) « 27{p - p c ) d ^/ 2 (dashed line in Fig. |), 
suggesting that u> c (p) scales as w^(p), but with a different 
prefactor (ui w 27 instead of wq ~ 6). 
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FIG. 3.: Phase diagram of vibrational modes of percolation 
clusters in d — 3. The FPC at UJ^(p) is marked by the contin- 
uous black line. The black discs indicate the LDT at ui c (p)- 
Finite size effects do only occur in the shaded region at very 
low frequencies. While the vibrational modes are (extended) 
phonons for low frequencies ui < u>((p) and localized frac- 
tons for uj > w c (p), a new intermediate regime of extended 
fractons is located between the two curves. 
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Next we consider d = 2. The results of the level statis- 
tics for vibrational modes on incipient infinite percolation 
clusters at p = 0.75 are shown in the inset of Fig. |[ We 
find no indication of a LDT. For u> > 0.15, the modes 
turn out to be localized since Io( w ) increases with in- 
creasing system size. This limit is well below the FPC 
frequency to^ ks 0.43 that we obtained from Eq. (Q) with 
too ~ 15. Hence, we can conclude that the vibrational 
modes remain localized even for frequencies well below 
the FPC frequency wg. For very small frequencies to, the 
level statistics does not provide a conclusive picture due 
to the restricted system sizes. Our results for the local- 
ization behavior in d = 2 are summarized in the main 
part of Fig. ^. 
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FIG. 4.: Phase diagram of vibrational modes of percola- 
tion clusters in d — 2. The black line separates localized 



fractons with ~ uo 



-2/d, 



from localized phonons with 



~ uj~ . In the shaded region, the method is not spe- 
cific since finite size effects occur. The inset shows 7o(w) for 
site concentration p — 0.75. The symbols correspond to dif- 
ferent system sizes, L = 70 (□), L = 100 (O), L = 150 (A), 
L = 220 (y), L = 350 (O), and L = 500 (+), correspond- 
ing to quadratic systems with linear sizes ranging from 95 to 
680. There is no indication of extended states. 

In conclusion, we have applied level statistics to inves- 
tigate the localization behavior of vibrational excitations 
on the infinite percolation cluster in d — 2 and d = 3. 
We have found that above p c two length scales have to 
be distinguished, (i) The wavelength A is strongly re- 
lated to the density of states g(u>) and its scaling can 
be derived from the analogy to diffusion. At wg, both 
X(u>) and g{u>) cross over from fractal to normal behav- 
ior, (ii) The localization length is supposed to diverge 
at the localization-delocalization transition, which occurs 
in d = 3 at a critical frequency w c , that is considerably 
larger than the fracton-phonon crossover frequency lo^. 
As a consequence, a novel regime of extended fractons 
exists in the phase diagram of vibrations in d = 3. These 
extended states might be observable in real experiments. 
In d = 2 there is no localization-delocalization transition, 



so also the phonon modes are localized. Thus, the vibra- 
tional modes in percolation systems show features char- 
acteristic of both classical diffusive and electronic sys- 
tems. For a complete understanding of the problem both 
analogies are useful and necessary. 

We thank the Deutsche Forschungsgemeinschaft for fi- 
nancial support. 



[1 
[2 

[s: 

[e: 
[7: 
[a: 

[9 
[10 

[11 
[12: 
[13; 

[14 

[15: 

[ie; 

[17 

[is; 

[19 

[20" 
[21 

[22" 
[23; 

[24; 

[25; 

[26; 



D. Stauffer and A. Aharony, Introduction to Percolation 
Theory, 2nd ed. (Taylor & Francis, London, 1992). 

A. Bunde and S. Havlin (eds.), Fractals and Disordered 
Systems, 2nd ed. (Springer Verlag, Heidelberg 1996). 

S. Alexander and R. Orbach, J. Phys. (Paris) Lett. 43, 
L-625 (1982). 

R. Orbach, Science 231, 814 (1986). 

S. Havlin and D. Ben-Avraham, Adv. Phys. 36, 695 

(1987). 

G. S. Grest and I. Webman, J. Phys. (Paris) Lett. 45, 
L-1155 (1984). 

B. W. Southern and A. R. Douchant, Phys. Rev. Lett. 
55, 966 (1985). 

K. Yakubo and T. Nakayama, Phys. Rev. B 40, 517 

(1989) . 

P. de Vries, H. de Raedt, and A. Lagendijk, Phys. Rev. 
Lett. 62, 2515 (1989); H. E. Roman, S. Russ, and A. 
Bunde, ibid 66, 1643 (1991). 

C. J. Lambert and G. D. Hughes, Phys. Rev. Lett. 66, 
1074 (1991). 

A. Bunde, H. E. Roman, S. Russ, A. Aharony, and A. B. 

Harris, Phys. Rev. Lett. 69, 3189 (1992). 

A. Petri and L. Pietronero, Phys. Rev. B 45, 12864 

(1992); A. Petri et al, ibid 49, 15067 (1994). 

P. Argyrakis, S. N. Evangelou, and M. Magoutis, Z. Phys. 

B 87, 257 (1992); S. N. Evangelou and P. Argyrakis, 

Phys. Rev. B 51, 3489 (1994). 

J. W. Kantelhardt and A. Bunde, Phys. Rev. E 56, 6693 
(1997). 

Y. Tsujimi, E. Courtens, J. Pelous, and R. Vacher, Phys. 
Rev. Lett. 60, 2757 (1988); R. Vacher et al, ibid 6b, 1008 

(1990) . 

M. Montagna et al, Phys. Rev. Lett. 65, 1136 (1990). 

E. Stoll, M. Kolb, and E. Courtens, Phys. Rev. Lett. 68, 
2472 (1992). 

S. Russ and B. Sapoval, Phys. Rev. Lett. 73, 1570 (1994). 
T. Nakayama, K. Yakubo, and R. L. Orbach, Rev. Mod. 
Phys. 66, 381 (1994). 



P. Grassberger, preprint bond-mat/9808095| (1998) 

R. Rammal and G. Toulouse, J. Phys. (Paris) Lett. 44, 

L-13 (1983). 

B. I. Shklovskii et al, Phys. Rev. B 47, 11487 (1993). 
E. Hofstetter and M. Schreiber, Phys. Rev. Lett. 73, 3137 
(1994). 

M. Batsch, L. Schweitzer, I. Kh. Zharekeshev, and B. 

Kramer, Phys. Rev. Lett. 77, 1552 (1996). 

M. L. Mehta, Random Matrices, 2nd ed. (Academic 

Press, San Diego, 1991). 

K. B. Efetov, Adv. Phys. 32, 53 (1983). 



4 



